Distribution and properties of fragments and debris from the split 
comet 73P/Schwassmann-Wachmann 3 as revealed by Spitzer 

Space Telescope 

William T. Reach 

Infrared Processing and Analysis Center, MS 220-6, California Institute of Technology, 

Pasadena, CA 91125 

reach@ipac . caltech . edu 

Jeremie Vaubaillon 

Infrared Processing and Analysis Center, MS 220-6, California Institute of Technology, 

Pasadena, CA 91125 

Michael S. Kelley 

Department of Astronomy, University of Minnesota, Minneapolis, MN 55455 

Carey M. Lisse 

Johns Hopkins Applied Physics Laboratory, SD/SRE-MP3/E167, 11100 Johns Hopkins 

Road, Laurel, MD 20723 

Mark V. Sykes 

Planetary Science Institute, 1700 E Ft. Lowell, Suite 106, Planetary Science Institute, 

Tucson, AZ 85719 

Received ; accepted 



Manuscript Pages: [79], Tables: EJ Figures: [15] 



- 2 - 

Proposed running head: Debris and Fragments of Comet 73P 

Editorial correspondence to: 

Dr. William T. Reach 

Infrared Processing and Analysis Center 

MS 220-6 

Caltech 

Pasadena, CA 91125 
Phone: 626-395-8565 
Fax: 626-432-7484 
E-mail: reachipac.caltech.edu 



- 3- 



During 2006 Mar - 2007 Jan, we used the IRAC and MIPS instruments on the Spitzer 
Space Telescope to study the infrared emission from the ensemble of fragments, meteoroids, 
and dust tails in the more than 3 degree wide 73P/Schwassmann-Wachmann 3 debris field. 
We also investigated contemporaneous ground based and HST observations. In 2006 May, 
55 fragments were detected in the Spitzer image. The wide spread of fragments along the 
comet's orbit indicates they were formed from the 1995 splitting event. While the number 
of major fragments in the Spitzer image is similar to that seen from the ground by optical 
observers, the correspondence between the fragments with optical astrometry and those seen 
in the Spitzer images cannot be readily established, due either to strong non-gravitational 
terms, astrometric uncertainties, or transience of the fragments outgassing. The Spitzer 
data resolve the structure of the dust comae at a resolution of ~ 1000 km, and they reveal 
the infrared emission due to large (mm to cm size) particles in a continuous dust trail 
that closely follows the projected orbit. We detect fluorescence from outflowing CO2 gas 
from the largest fragments (B and C), and we measure the CC^H^O proportion (1:10 and 
1:20, respectively). We use three dimensionless parameters to explain dynamics of the solid 
particles: the rocket parameter a is the reaction force from day-side sublimation divided by 
solar gravity, the radiation pressure parameter (3 is the force due to solar radiation pressure 
divided by solar gravity, and the ejection velocity parameter v is the particle ejection speed 
divided by the orbital speed of the comet at the time of ejection. The major fragments have 
v > a > (3 and are dominated by the kinetic energy imparted to them by the fragmentation 
process. The small, ephemeral fragments seen by HST in the tails of the major fragments 
have a > v > (3 and are dominated by rocket forces (until they become devolatilized). The 
meteoroids along the projected orbit seen by Spitzer have (3 ~ v ^> a and are dominated by 
radiation pressure and ejection velocity, though both influences are much less than gravity. 
Dust in the fragments' tails has (3 ^> [y + a) and is dominated by radiation pressure. 

Keywords: comets, meteors, infrared observations 
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Introduction 



Comet 73 P / Schwass mann- Wachman n 3 split in 1995, with major fragmentation 



beginning on 1995 Sep 7 (ISekanina 



20051 ) . Phenomena associat ed with the split included 



dramatically increased OH gas emission (ICrovisier et al 



1996), an optical brightness 



increase by 4 magnitudes, and the presence of at least two and likely 4 fragments. Two 
of the fragments (73P-B and 73P-C) from the 1995 event were long-lived and returned on 
both the 2001 and 2006 perihelion passages. (Based on the comet's 5.39 yr period, the 
next perihelion passage will be in 2011 Oct.) 73P-C is the brightest fragment and 'leads' 
73P-B in the orbit, suggesting it is the principal and most massive remnant of the pre-1995 



comet, 



jy comparison to other split comets where the main fragment generally leads the 



others ( Boehnhardt 



20041 ). In the 1990 apparition, prior to the split, 73P was observed as 
a single, intact comet, with an upper limit on fragments being 3.2 mag fainter than the 



primary (IChen and Jewitt 



19941 ) . Yoshida's compilation of the lightcurvq^l shows a typical 
brightness profile for a Jupiter-family comet in 1990 (with peak heliocentric magnitude 
11.5), and a nearly identical brightness trend b eing followe d in 1 995 up until 1995 Aug, 



Lisse et al. 



(119981 ) reported significantly 



when it brightened by 6 magnitudes. However, 
enhanced mid-infrared emission from the comet detected by the Diffuse Infrared Background 
Experiment on the Cosmic Background Explorer in 1990, suggesting increased emission of 
heavy dust particles. On the 2001 return, fragments B and C were observed and followed a 
typical activity curve, elevated by approximately 3.5 mag with respect to that observed in 
1990. At least one new fragment, 73P-E, was found in the 2001 return but it quickly faded 
from view (and did not return in 2006). 

In 2006, as the fragments of 73P again approached perihelion, there was renewed 
interest due to 73P's anticipated close approach to the Earth (A = 0.08 AU) on 2006 May 



x http: / / www.aerith.net / comet / catalog/0073P 
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13, just before perihelion (q = 0.939 AU) on 2006 Jun 7. That the comet's orbit nearly 
intersects the Earth's is significant not only because it becomes bright for ground-based 
observations but also because the comet has an associated meteoroid stream, the r Herculids 



( Jenniskens 



20061 ). and the comet is also an excellent mission target, as it can be reached 
relatively inexpensively by spacecraft. (Indeed 73P was considered as a target for the 
ROSETTA and CONTOUR missions.) Amateur and professional observations revealed an 
extensive set of fragments along the orbit of 73P, with the brightest fragment still being 
73P-C (though in late 2006 Apr, 73P-B was brighter) and some small fragments found 
ahead of 73P-C. Yoshida's compilation of lightcurves shows a typical shape for 73P-C, 
though it remained some 2.5 mag brighter than 73P was in 1990. The lightcurve of fragment 
B has at least one significant increase, becoming some 2.5 mag brighter beginning around 
2006 Mar 31 and returning to its pre-flare-up trend by 2006 May 8. Fragment G was first 
observed in 2006 Feb and it had a significant flare-up by approximately 2 magnitudes, 
beginning 2006 Apr 5 then fading from view in early May. The flare-ups of fragments B 



and G were almost certainly ass ociated wit 



by ground-based te 
HST and Chandra; 



escopes (e.g. 



Weaver et al. 



i sign ificant sub-fragmentation as was observed 



use et al. 



2006 . 



2007) and sp ace-based observatories (including 



Wolk et al. 



20091 ). 



By studying the properties of the fragments of split comets, we hope to determine 
the properties of fresh surfaces of comet nuclei as well as the structure of the nucleus. 
Comet nuclei are heated to depths of several cm during diurnal (and of order meters for 
orbit-averaged) heat pulses, and they lose of order 1 m of their surface each perihelion 
passage. When comets split into large fragments, fresh interior surfaces are exposed, 
allowing observation of relatively primitive material from the epoch of the comet's formation. 



The major fragments B and C were observed spec troscopically anc 



to each other, both being carbon-chain depleted (jKobayashi et al. 



found to be very similar 



20071 ). The similarity 



of their composition is in contrast to the fact that fragment B was actively producing 
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icy sub-fragments during the 2006 apparition while fragment C appeared to exhibit only 
'normal' activity without sub-fragmentation in 2006. Thus the chemical properties of the 
ice may not be very different at depth than they are for active regions in the surface, 
but the distribution of ice within the nucleus may be highly nonuniform. In addition 
to the chemistry of the freshly-outgassing surfaces, we can learn about the structure of 
cometary nuclei from the way they split. Fragmentation will proceed rapidly for sub units 
that are weakly bound, possibly leading down to the size of cometesimal building blocks. 
Debris of a wide range of sizes was inferred from observations of C/1999 S4 (LINEAR), 
with comparable mass in large (~ 0.3 m) meteor oids and detectable (~ 40 m) fragments 



( jWeaver et al. 



2001 



Bockelee-Morvan et al. 



200lh . 



In this paper we present Spitzer Space Telescope mid-infrared observations of the 
fragments and debris from 73P. The outstanding sensitivity of space-based mid-infrared 
observations for extended, inner solar system objects makes these observations unique in 
their ability to simultaneously image the fragments, meteoroids, and dust. In a companion 
paper, we analyze the dynamics of the dust and meteoroids to infer the size distribution, 
finding that there is an excess in the abundan ce of mm-sized and larger pa rticles relative 



to a power-law fitted to the smaller particles (jVaubaillon and Reach 



20081 ). This paper 



concentrates on the physical properties and dynamics of the fragments. 



2. Observations 

More than 57 detected discrete nuclear fragments, associated meteoroids, and the 
anti-solar dust tails of 73P/SW-3 span a complex debris field along more than 3 degrees of 
mean anomaly of the comet's orbit. This debris field was mapped three times as part of 3 
Spitzer observing programs. All observations were designed to cover the region spanning 
from ahead of fragment C to well behind fragment B. Table [1] lists the observing conditions, 
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including the length of the arc along the orbit covered by each mosaic. To illustrate the 
viewing geometry at the times of the Spitzer observations, Figure [T] shows the locations of 
Earth, Spitzer, and the comet. 



2.1. 2006 Mar-Apr 



The Spitzer Cycle 2 GO imaging observations (Pro gram ID 20039) in 2006 Mar-Apr 



were performed with the Infrared Array Cam era (1RAC 



band Infrared Photometer for Spitzer (MIPS 



Rieke et al. 



azio et al 



2004 ) and the Multi- 



20041 ). with very good viewing 



conditions — the comet was close to Earth and Spitzer, and the line of sight was inclined by 
20° from the orbit plane (in contrast to common comet observing conditions at orbit plane 
angles < 5°). The map center and size were tuned based on rapidly evolving knowledge of 
the fragment locations as of 2006 February, as gleaned from amateur astronomers' internet 
postings and our own reconnaissance using the Palomar Observatory Hale telescope. The 
images were made twice, with a separation of 5.4 hr between epochs; during the time 
between epochs, the comet moved by 14.9', so a completely different background field of 
galaxies was present. 

We combined the images at the two epochs into a comet-only image as follows. First, 
we registered the mosaics made for each epoch, so that the comet fragments appear in 
the same pixels. Then, we created a celestial-background-subtracted image by taking, for 
each pixel, the lowest of the two observed brightnesses (epochs 1 and 2). This method 
works because the sky brightness is the sum of comet emission plus only positive emission 
from galaxies and stars. The minimum between two images shifted by 14.9' (in celestial 
coordinates) eliminates all non-cometary emission on angular scales smaller than 14.9'. 
There is some diffuse, extended interstellar medium emission in the field, which is not 
completely removed, but its brightness is much smaller than the comet emission. Also, there 



- 8- 



Table 1. Log of observations 



Start Date 


R 


A s a 


Wavelengths 


Arc b 


Tiles 


Ttile C 


5T shad A 


N- ° 

iv image 


(UT) 


(AU) 


(AU) 


(Mm) 


(10 6 km) 


(N x ' x ') 


(min) 


(min) 




2006-03-28 09:35 


1.37 


0.64 


3.6,4.5,5.8,8 


3.4 


3 x 47 x 30 


40 


122 


1,854 


2006-04-01 19:19 


1.33 


0.60 


24, 70, 160 


3.6 


1 x 124 x 82 


181 


325 


5,200 


2006-05-04 12:49 


1.06 


0.37 


24, 70, 160 


6.7 


4 x 131 x 84 


181 


181 


20,800 


2007-01-01 17:16 


2.80 


2.58 


24 


5.2 


1 x 46 x 8 


162 




672 



a Distance from the observatory (Spitzer) to the comet 
b Distance along the comet's orbit covered by the observation 
c Duration of one map tile 

d Interval between primary and shadow observations (for each tile) 

c Number of images (5' x 5' field of view for IRAC at each wavelength or MIPS at 24 /zm) 
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Fig. 1. — The relative positions of Spitzer (□), Earth (©), and comet 73P (A) on the 
dates of our MIPS observations. The positions are shown in a geocentric ecliptic coordinate 
system, and the orbits of Spitzer and 73P are plotted as gray lines. 
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are rare cases where, by chance, galaxies overlap between the images at the two epochs; we 
verified none of the comet fragments are galaxies by inspecting the single-epoch images for 
each fragment. A minor caveat to the method is that noise is both positive and negative 
(relative to the mean), so our method will generate an image with a systematic negative 
bias. The noise is primarily due to shot noise from the sky brightness, so the negative 
bias is proportional to the square-root of the sky brightness. The signal-to-noise of the 
images is high (> 100) due to the bright zodiacal background, so the negative bias is small: 
< 1% estimated from numerical simulation. It is important to note that the background 
subtraction required no assumption about the cometary emission, so cometary emission on 
all angular scales was preserved. The most prominent artifacts in the comet-only image are 
5'-wide vertical stripes; these are caused to gradual detector drifts, and their shape on the 
image is due to the scanning strategy. 

Figure [2] shows the MIPS image. The debris trail is clearly evident, spanning the entire 
image from ahead of fragment C to well behind fragment B. The MIPS arrays were scanned 
across the field in 26 fast scan legs with full-array (5') spacing between legs. A total of 6 hr 
of observing time was spent on this image. 

Figure [3] shows the IRAC image. Because the fragment and debris field of the comet 
spanned such a large field on the sky (exceeding the maximum observation size), the image 
was created by observing three tiles staggered along the orbit of the comet. For each tile, a 
raster map of 10 x 6 pointings spaced by one 5' x 5' IRAC field of view was made, and two 
12 sec frames were taken at each pointing. A total of 4.4 hr observing time was spent on 
this image. 
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Fig. 2. — MIPS 24 /j,m image from 2006 Apr 1. Major fragments C, B, and G are labeled, 
as are the smaller fragments A2, A3, and A4. North is up, East is to the left, and the image 
size is 123' x 83'. The image was background-subtracted, using the minimum between the 
two observing epochs as described in the text. 



Fig. 3. — IRAC color-composite image from 2006 Mar 28. Blue is 3.6 /im, green is 4.5 /im, 
orange is 5.8 /im, and red is 8 /im. North is up, East is to the left, and the image size is 
110' x 63'. The diffuse green emission, extending sunward from the main fragments, is due 
to CO and C0 2 gas emission (JTJ). The image was are background-subtracted, using the 
minimum between the two observing epochs as described in the text. 
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2.2. 2006 May 

The Director's Discretionary observations (Program ID 274) in 2006 May were 
performed with MIPS and spanned a larger area both on the sky and in the comet's frame, 
with excellent viewing conditions — the comet was at its closest approach to Spitzer and the 
orbit plane angle was 25°. The observations were designed and proposed immediately after 
receipt of the 2006 Apr images were received. The field size for the 2006 May observations 
was calculated so as to cover the orbit from well ahead of fragment C to behind the 
by-then-discovered fragment G (which was near the edge of the 2006 Apr observation), 
after adjusting for the relative distance between the comet and the observatory on the two 
dates. It is worth noting that the viewing geometry from Spitzer was somewhat different 
from that for ground-based observations. As shown in Figure [lj the comet was much closer 
to Earth (0.10 AU) than it was to Spitzer (0.37 AU); in fact, the comet was closer to Earth 
than Spitzer was to Earth (0.32 AU). 

Figure H] shows the MIPS image. To cover the large field (larger than that which can 
be covered in a single scan map) the field was broken into 4 rectangular tiles, with centers 
aligned along the projected orbit including overlaps at the corners. Each tile is the same 
type of observation as the 2006 Apr image. A total of 24 hr of observing time was spent on 
this image. 

2.3. 2007 Jan 

The Spitzer Cycle 3 Guaranteed Time observations (Program ID 30010) in 2007 Jan 
were designed to cover a similar portion of the orbit, but at much less favorable viewing 
geometry. As shown in Figured], the comet was much more distant; further, the orbit plane 
angle was only 2.3°, making the trajectories of small and large particles overlap on the sky. 
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Fig. 4. — MIPS image from 2006 May 4-6. These images are background-subtracted, using 
the minimum between the two observing epochs as described in the text. The orientation is 
indicated with arrows pointing N and E, and the length of those arrows is 17'. Some major 
fragments are labeled. Panel (a) is the is the entire mosaic, combining the 4 large tiles. 
Panels (b) zooms into a region behind 73P-G that is particularly rich in fragments. Panel 
(c) zooms into the portion of the observation furthest behind the main fragments, and it 
includes the fragment located furthest perpendicular to the project orbit. 
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The observation was made using the MIPS 24 /im photometry mode, with one cycle of 
10-sec-exposure photometry performed at each of a set of cluster positions spaced along the 
projected orbit. 

Figure [5] shows the MIPS image. The anti-solar and projected comet motion vectors 
are nearly aligned on the sky, so the tail (small particles fro m the current orbi t) and trail or 



Reach et al 



20071 ) overlap the 



meteoroid stream (larger particles from previous orbits, cf. 
fragment distribution. Though the heliocentric distance was 2.8 AU, and the sublimation 
rate of H2O ice should be relatively low, the large and small particles cannot be readily 
disentangled, as they could in the 2006 Apr-May observations. 

Only fragments C and B are detected, despite having enough sensitivity to have 
detected fragment G if it scaled similarly to the other main fragments. The lack of a 
detectable condensation at the predicted location of fragment G suggests this fragment 
completely disrupted. This is consistent with the appearance of fragment G in Hubble Space 
Telescope images: instead of a dominant nucleus with retinue of smaller bodies, fragment 
G appeared as a swarm of smaller bodies. 

The diffuse emission in 2007 Jan closely follows the projected orbit. A debris trail 
stretches from the one edge of the image to the other, i.e. from well ahead of fragment C 
to well behind fragment B. The emission is far brighter than a typical debris trail at 2.8 
AU from the Sun. Ahead of fragment C, the diffuse emission is at 0.1 MJy sr _1 , which 
is typical of debris trails produced natu rally by comets a nd representing mass loss on 



previous revolutions around the Sun (cf. 



Reach et al. 



20071 ). But following fragments C and 



especially B, the brightness is much higher and has an atypical brightness distribution, as 
we will discuss further in £18.21 below. 
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Fig. 5. — MIPS image from 2007 Jan. The image is rotated for compact display; the 
orientation of celestial North and East are indicated. The image covers 46' x 18'. Fragments 
C and B are labeled, as is the expected location of fragment G if it had remained at the 
same spacing behind fragment B, relative to the separation between fragments B and C, in 
2007 Jan as it had in 2006 Apr and 2006 May. 
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3. Properties of fragments 

3.1. Spectral energy distributions of the comae of fragments B and C 

Figure [6] and Table [2] show the spectral energy distribution of fragments B and C. 
Fluxes were measured within a 25" radius aperture, and they are dominated by dust in 
the coma. Fragment B was measured at all wavelengths (3.6-70 /mi), but fragment C was 
saturated at 24 /tm and was not covered by the operative portion of the 70 /im detectors. 
To estimate the 24 /mi flux of fragment C, we replaced the central 3-pixel (7.5") radius 
region centered on the nucleus (the only pixels that are saturated) with a 1/p brightness 
distribution, convolved with a 5" gaussian FWHM beam and fitted to the coma from 
10-30" from the nucleus. The estimated flux in the saturated portion of the image is a 
small fraction (< 10%) of the total; nonetheless, we adopt a conservative 40% uncertainty. 
Note that the observations at wavelengths shorter than 10 /mi were taken on 2006 Mar 
28, while those at longer wavelengths were taken on 2006 Apr 1. Fragments B and C 
were brightening through 2006 Mar and Apr. We fit greybody, F v oc B V (T), and modified 
greybody F v oc u a B u (T), where T is the temperature and B v is the Planck function, to the 
spectral energy distributions, allowing for a fractional brightening over the 4-day interval, 
between 2006 Mar 28 and 2006 Apr 1, in the range 1-1.4. If the brightness increase is more 
than 1.4, then it is not possible to put a modified blackbody through the 5.8-24 /tm data 
points that are completely dominated by thermal dust emission. 

The spectral energy distribution of fragment B could not be adequately fit using 
a greybody spectrum, irregardless of any possible brightness variation between the two 
observing dates. The observed 70/24 /tm flux ratio of 0.13 ± 0.02, on 2006 Apr 1, requires a 
color temperature > 800 K, which is far higher than comet grains at 1.33 AU are expected 
to be, and which is completely inconsistent with the color temperature from of 255 ± 12 
K measured from the 5.8-8 /tm brightness on 2006 Mar 28. A modified blackbody yields 
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an acceptable fit to the spectrum, as long as the emissivity of the grains decreases as the 
—0.79 ± 0.12 power of the wavelength. The grain temperature from such a fit is 236 ±8 K, 
which is consistent with the 238 K expected for rapidly-rotating isothermal grains at 1.37 
AU from the Sun. 

The depressed 70 /zm brightness is likely due a dominant contribution from grains too 
small to emit efficiently at such long wavelengths. Particles smaller than \/2ir = 11 /im are 
expected to be in the small particle limit, where their emissivity will decrease as A^ 1 . The 
observed emissivity index of 0.79 suggests the dominant particles are in the range 2 to 20 
/iin in radius. 
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10 

Wavelength (/xm) 



Fig. 6. — Spectral energy distribution of fragments B and C based on observations on 
2006 Mar 28 and Apr 1. The dashed lines for both fragments show a fit with a = and 
free temperature, fitted only to the 5.8-24 /iin photometry. These fits clearly overpredict 
the long-wavelength flux. For fragment B the thick solid line shows a best fit with a free 
temperature and emissivity index (see text for definition), with best-fit values a = 0.79±0.12 
and T = 236 ±8 K. For fragment C the thick solid line shows the best fit with fixed a the 
same as found for fragment B, and free temperature T = 239 ± 4. 
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Table 2. Fluxes of fragments B and C on 2006 Mar 28 a 



Wavelength 


Fragment C Flux (Jy) 


Fragment B Flux (Jy) 


5.65 


0.37 ±0.02 


0.115 ±0.01 


7.74 


2.1 ±0.1 


0.63 ±0.05 


23.7 


9.6 ±2 


3.0 ±0.2 


71.4 




0.40 ±0.03 



a Fluxes are measured within a circular 25" radius aperture. The 23.7 and 81.4 /im fluxes 
were scaled from 2006 Apr 1 measurements to match the 2006 Mar 28 observations at the 
other wavelengths. 
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3.2. Locations and fluxes of intermediate-sized fragments 

Table [3] lists the fragments detected in the 24 /zm images at each of the three observing 
epochs. Only the primary fragments B and C, which were most likely produced in the main 
1995 splitting event, are present in all three epochs. Fragment G was very bright before 
the 2006 perihelion passage but was not recovered afterward. Our 2007 Jan observation 
shows no nucleus or coma at the location of fragment G. The 2006 May observation shows 
numerous fragments, as did ground-based optical observations around the same time. 
As discussed later in this paper, the infrared and optical fragments could not be readily 
matched. We used the JPL Horizons utility to calculate locations of all fragments for 
which orbits had been calculated (based on ground-based optical astrometry). Fragment G 
corresponds well with the third-brightest mid-infrared fragment (2006 Apr and May), and 
the optically-bright fragment R's predicted location is close to an infrared-bright fragment 
in 2007 May. These identifications are noted in Table [31 Three other fragments whose 
locations in the Spitzer image are close to predicted locations of fragments with calculated 
orbits are listed with a question mark in Table [31 The remaining fragments have no clear 
counterpart. We will refer to fragments in the 2006 Apr image as A1-A5, and those in 
the 2006 May image as S1-S55, except for the well-identified fragments. It is likely that 
fragment A4 and S16 (located just ahead of fragment B) are the same object. 

The flux of each fragment was measured using annular-background-subtracted aperture 
photometry. Fragments B and C are saturated in the images so their fluxes could not be 
measured. Table [3] lists the fluxes as well as a note on the presence or lack of a dust tail, by 
which we mean diffuse emission that extends from the fragment in the anti-solar direction. 
The observed flux of a fragment combines information about the nucleus size and the 
current dust production rate. For bare nuclei, the flux at 24 fim is FiD 2 where Fi is the flux 
of a body with a diameter of 1 km. We calculated F\ = 88 mJy for the asteroid standard 
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thermal model (ILebofsky et a/.lll986l ). 90 mJy for a rapidly-ro tating (isothermal 



model, and 42 mJy for the near-Earth-asteroid thermal model ([Harris et al 



-latitude) 



19981 ). assuming 



albedo 4% and the 2006 May observing geometry from Spitzer. For these calculations it is 
important to note the high phase angle (Sun- Comet- Spitzer) of 72° is outside the normal 
range for main belt asteroids. We adopt F\ = 50 mJy for diameter estimates this paper. 
Assuming an albedo of 0.04, the size of fragment C was measured using Hubble Space 



Telescope to be 0.68 ± 0.04 km flToth et al 



20051 ); the nuclear flux at 24 /im would then 



be 23 mJy. The observed 24 /im flux is 2 orders of magnitude higher than this value, 
suggesting it is dominated by thermal emission from dust surrounding the nucleus. 

Many of the smaller fragments have tails, showing they, too, are actively producing 
dust. The observed fluxes yield only upper limits to the fragment sizes; the upper limits 
can be obtained from the fluxes, F, in Table [3] as D < (F/50 mJy) - 5 km. The faintest 
fragments, with fluxes ~ 3 mJy have a nuclear diameter limits < 250 m. These limits are 
likely much larger than the actual fragment sizes, since the fragments are active and their 
mid- infrared fluxes are much smaller than that of fragments C and B. 

The distribution of observed fluxes gives some indication of the size distribution, if the 
dust production rate scales with nuclear size. In fact the dust production is not expected 
to be perfectly correlated with size and may even be episodic, so keep in mind that the 
following analysis of fragment 'size' distribution may be skewed relative to the nucleus size 
distribution. However, the fluxes (from aperture photometry) measured in the infrared 
images should be directly proportional to those measured from optical images, so that 
we can compare the Spitzer results to previous optical measurements. Figure [7] shows 
the cumulative luminosity distribution of the 53 fragments with fluxes measured 2006 
May. We fit pieces of the cumulative luminosity distribution with a power-law function to 
approximate the true cumulative luminosity function. Our results are presented in Table HI 
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The cumulative luminosity function has a clear break in the power-law slope at a flux of 10 
mJy. 

The s 



Fuse et al. 



ope fo r fragments fainter than 10 mJy (-0.42) is close to the slope found by 
(120071 ) (-0.555) in their Subaru observations of mini-fragments near 73P-B. The 
mini-fragments near 73P-B fall along its tail and are likely icy (as discussed in §5]); the 
Spitzer-detected fragments fall along the comet's orbit and likely have less exposed ice, as 
explained in detail in $6] below. The Spitzer-detected fragments are also likely somewhat 
larger (~ 50 m) than the Subaru-detected mini- fragments (~ 20 m). The similarily of 
flux distribution slopes for the fragments along the trail (from the breakup in 1995) and 
the mini-fragments in 73P-B's tail (from the breakup of 73P-B in 2006) suggests that the 
fragmentation process itself operates similarly over a wide range of size scales. 

If we interpret the luminosity distribution slopes as nuclear size distributions, where 
the size distribution is 

n(R)dr oc R~ q dr, (1) 

then the inferred size distributions slopes are q = 1.84 for small (F < 10 mJy) fragments 
and q = 2.56 for large (F > 10 mJy) fragments. We compare our derived size distributions 
to other comet size distributions in Table [51 For a system in collisional equilibrium, q = 3.5, 
which is much steeper (i.e. more, smaller particles) than we observe and steeper than 
any of the comet size distributions in Table O Our large fragment size distribution is 
similar to that derived from three different studies: (1) the size distribution of fragments 
discovered after the c atastrophic breakup of comet C/1999 S4 (LINEAR), q = 2.72 ± 0.60 



(IMakinen et al 



20011 ) . (2) the inferred in itial size distribution of Kreutz group sungrazing 



comets (R = 20 — 200 m; ISekanina 



comet nuclei for R > 1.6 km (ILamy et al. 



2003), and (3) the observed size distribution of ecliptic 



20041). 



These agreements suggest that the break in the slope at F 2 4 = 10 mJy indicates a 
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change in the size distribution slope, rather than a change in the dependence of coma 
brightness on fragment size. The change in size distribution slope doe snot necessarily 
reflect the fragmentation process; instead, the relative deficit of smaller fragments may 
indicate that they have a short lifetime. Fragments with a short lifetime will be more 
abundant in images taken soon after an event that produces small fragments. They will 
also be more abundant in higher-angular-resolution images taken close to the site of their 
production. Both of these conditions are met for the Subaru and HST observations, which 
detect more of the very small fragments than we could with Spitzer. 
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Table 3. Fluxes and locations of Fragments of 73P/Schwassmann-Wachmann 3 



Fragment F 2A D xC /D B c & D xC /A h z/A c Comment 
(mJy) (') (') 



2006 Apr 



c 


4760 











tail> 70' 


A2 


4.5 


-0.63 


-61.3 


< 0.1 




A3 


1.9 


-0.81 


-79.1 


< 0.1 


fuzz 


A4 


10.8 


-0.95 


-93.3 


< 0.1 


fuzz 


B 


2690 


-1 


-98.2 


< 0.1 


tail> 18' 


G 


138 


-1.18 


■115.8 


< 0.1 


coma+tail 


2006 May 


SI 


2.3 


1.11 


149.1 


2.23 


fuzz 


S2 


4.7 


0.96 


129.1 


1.93 


fuzz 


BR? 


58.9 


0.79 


105.7 


1.25 


tail 3' 


AT? 


4.3 


0.71 


95.1 


0.79 




S5 


12.9 


0.60 


81.0 


0.76 


tail 3' 


S6 


5.2 


0.48 


64.8 


0.24 




AV? 


22.7 


0.47 


62.8 


0.08 


tail 3' 


S8 


3.6 


0.45 


61.1 


0.31 
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Table 3 — Continued 



Fragment F 2A D xC /D B c & D xC /A h z/A c Comment 
(mJy) (') (') 



S9 


70.1 


0.27 


36.8 


0.32 


tail 4' 


S10 


25.2 


0.26 


35.3 


0.25 


tail 3' 


c 




0.00 


0.0 


0.00 


saturated 


S12 


38.1 


-0.55 


-74.0 


-0.40 


tail 2' 


S13 


13.8 


-0.60 


-80.8 


0.24 


tail 2' 


S14 


10.5 


-0.62 


-82.7 


-0.08 


tail 2' 


S15 


4.7 


-0.86 


-116.0 


0.16 




S16 


3.2 


-0.95 


-127.7 


-0.08 


A4? tail 1.5' 


B 




-1.00 


-134.3 


0.28 


saturated; tail 90' 


G 


336.9 


-1.19 


-159.5 


0.21 


tail 30' 


S19 


13.7 


-1.22 


-163.6 


0.59 


fuzz 


R 


219.3 


-1.29 


-172.7 


-0.10 


tail 5' 


S21 


47.8 


-1.30 


-174.7 


0.57 


BN? coma 


S22 


2.9 


-1.33 


-178.2 


-0.00 


fuzz 


S23 


2.9 


-1.34 


-179.8 


-0.13 


fuzz 


S24 


5.7 


-1.36 


-183.3 


0.40 




S25 


5.7 


-1.37 


-183.5 


0.15 




S26 


3.3 


-1.37 


-184.2 


0.95 




S27 


55.5 


-1.38 


-185.1 


0.49 


N??; tail 2' 


H 


64.4 


-1.39 


-186.3 


-0.05 


tail 2' 
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Table 3 — Continued 



Fragment F 2A D xC /D B c & D xC /A h z/A c Comment 





(mJy) 




(') 







S29 


13.7 


-1.39 


-187.4 


0.67 


BQ?? 


S30 


3.7 


-1.40 


-187.7 


-0.03 




S31 


27.6 


-1.41 


-188.9 


0.15 


BP?? 


S32 


9.9 


-1.43 


-192.7 


-0.01 




S33 


10.5 


-1.44 


-193.1 


0.16 




S34 


15.7 


-1.45 


-194.6 


0.33 




S35 


6.9 


-1.48 


-198.9 


0.67 




S36 


6.4 


-1.49 


-199.8 


0.13 




S37 


6.8 


-1.54 


-206.6 


0.25 




S38 


3.5 


-1.59 


-213.6 


0.02 




S39 


11.4 


-1.67 


-223.9 


0.38 


compact 


S40 


10.7 


-1.70 


-228.6 


0.29 




S41 


4.6 


-1.72 


-231.5 


0.54 




S42 


34.2 


-1.74 


-234.0 


0.82 


tail V 


S43 


4.8 


-1.76 


-236.7 


0.11 




S44 


19.1 


-1.77 


-237.9 


-0.38 


tail V 


S45 


11.0 


-1.78 


-239.3 


-0.53 




S46 


3.2 


-1.80 


-241.5 


0.04 




S47 


28.0 


-1.83 


-246.5 


4.16 


tail V 


S48 


26.4 


-1.84 


-246.9 


0.14 
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Table 3 — Continued 



Fragment F 24 D xC /D BC & D xC /A h z/A c Comment 
(mJy) (') (') 



S49 


43.7 


-1.85 


-247.9 


0.56 


S50 


12.4 


-1.86 


-250.1 


0.19 


S51 


8.8 


-1.88 


-252.3 


0.83 


S52 


7.5 


-1.91 


-256.3 


0.30 


S53 


7.3 


-1.96 


-263.8 


0.29 


S54 


14.9 


-2.05 


-275.5 


0.41 


S55 


10.8 


-2.05 


-275.8 


-0.08 



tail I' 



2007 Jan 

C 55.1 

B 13.2 -1 -18.1 

G <0.2 -1.18 d 

Reparation of fragment from 73P-C, in units of the separation between fragments B and 
Reparation of fragment from 73P-C, in arcmin 

c Perpendicular displacement of fragment from the projected orbital path of 73P-C, divided 
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by the Comet-Observer separation, in arcmin 

d assumed location of fragment, based on previous observations 
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Table 4. Best-fit power-law slope of the fragment cumulative luminosity function in May 

2006. 

Flux Range (mJy) k a a mad (dex) b 

F > 10 mJy -0.78 0.04 

F < 10 mJy -0.42 0.02 

a Power-law slope determined by the least absolute deviation method. 
b Mean-absolute deviation in log space. 



Table 5. Comet fragment differential size distributions. 



Comet 


Q 


Source 


73P (small) 


1.84 


this work 


73P (large) 


2.56 


this work 


73P-B tail fragments 


2.11 


Fuse et al. 2007 


C/1999 S4 (LINEAR) 


2.72±0.60 a 


Makinen et al. 2001 


Kreutz group 


2.7-3.0 b 


Sekanina 2003 


All comets 


2.9 ±0.3 


Lamv et al. 2004 



a Inferred from the fragment water production rates. 

b Derived from the cumulative mass distribution slope range, -0.56 to -0.68. 
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Fig. 7. — The cumulative 24 jum flux distribution of 56 fragments observed in May 2006. 
Three power-law fits to the luminosity di stribution are pre sented: the solid line shows the 



slope derived from Subaru observations by 



Fuse et al. 



(120071 ) ; the dotted line is a fit for fluxes 



less than 10 mJy; and the dashed line is a fit for fluxes larger than 10 mJy 
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4. Distribution of fragments 

Determining the identity and location of each of the fragments in our infrared image 
series is a non-trivial exercise. The location of each fragment was rotated into a coordinate 
system measuring separation from fragment C along the line directed toward fragment 
B, and the distance perpendicular to this line. To remove the geometric projection due 
to changing distance from the telescope, the distance along the orbit was divided by the 
separation between B and C. There are 10 fragments leading C, 5 fragments between B 
and C, and 38 fragments following B (for a total of 55 fragments). The infrared fragments 
span nearly the entire observed field, ranging from -2.07 to +1.13 times BC, where BC 
is the vector from B to C with origin at C. For comparison, orbits of 66 fragments have 
been generated by the Minor Planet Center and JPL/Solar System Dynamics Group. The 
similarity of the number of infrared and optical fragments suggests there may be good 
correspondence, and most of the infrared fragments may have optical counterparts. 

In an effort to cross-identify the fragments, we generated ephemerides using the 
Horizons service for all 66 fragments of 73P with known orbits. Figure [9] shows the 
distribution of fragments along the orbit of 73P in May 2006. There are 13 ephemeris 
positions ahead of C, 7 between C and B, and 44 behind B, The positions range from 
-2.95 to 1.05 times the separation between B and C, relative to C (except one fragment, 
S, predicted to be at 6.6 times further behind C than B). The distribution of predicted 
positions is very similar to the distribution of infrared fragments, especially after removing 
7 positions that fall outside the infrared image behind B. Good matches were found for the 
major fragments C, B, G, and R. Beyond that, to our surprise, few convincing matches were 
found. Likely matches include fragments BR and AT, for which the position coincidence is 
much better than expected by chance. Plausible matches include optically-bright M and N 
with moderately bright infrared fragments, though the deviations are large and multiple 
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matches are possible. Inspecting Figure [9] shows that part of the reason for the lack of 
matches is that the distance of the fragments perpendicular to the fragment C orbit is 
very different for the observed (Spitzer) and ephemeris positions. The root-mean-square 
perpendicular separation is 0.6' for the observed positions and 10' for the ephemeris 
positions. Typical quoted uncertainties in the predicted ephemeris positions are < 0.1'. 

Even if we ignore the out-of-plane motion, the locations of predicted positions along 
the orbit are different from those in the infrared image in detail. For each infrared fragment, 
we found the closest optical fragment along the orbit (neglecting perpendicular distance). 
The separation between the infrared and closest optical fragments has a distribution with 
a core dispersion of 0.005 |-BC|, containing ~ 30 fragments, with a wider distribution 
of width 0.03 \BC\. The core of this distribution could contain a significant number of 
matching objects, though its dispersion (40" on the sky) is broader than the resolution 
(5"). Fragments C, B, G, and R have infrared-optical matches within 5". Using a set of 
random fragment locations, the wider part of the separation distribution is reproduced. We 
have noted in Table [3] those fragments with optical matches within 5" as identifications in 
the first column; those within 10" as identifications with "?" in the first column; and those 
within 40" in the comment field with "?" . Figure [10] shows that the situation in 2006 Apr 1 
is similar to that in 2006 May 5: apart from the brightest fragments C, B, and G, there are 
no matches between the predicted positions and the other infrared fragments. 

An attempt to match the optical and infrared fragments, using a procedure very 
similar to the one described above, was undertaken by P. Birtwhistle (website and personal 
communication). One difference between our procedures was that, instead of using the 
Horizons-predicted ephemerides, Birtwhistle sorted fragments by time of perihelion and 
compared the Spitzer image to close-in-time optical images and astrometry reported to the 
Minor Planet Center. In general, finding matches was extremely difficult, as we reported 
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above. But there were some sections of the fragment train where matches appeared 
plausible. Figure [8] shows a dense section of the fragment train, with the optically-identified 
fragments labeled on the optical image and the infrared-identified fragments labeled on 
the infrared image. Both attempts at matching fragments are indicated. As noted above, 
the ephemeris positions do not often correspond with fragments in the Spitzer image. 
In Fig. [HI fragment 73P-H matches very well. The ephemeris positions of M and N are 
oriented similarly to two nearby fragments in the Spitzer image, and AN is close to a pair 
of fragments; however, none of these correspondences are convincing. The ordering and 
labeling by P. Birtwhistle allows possible identifications for M, BO, N, H, BQ, BP, and 
AN with Spitzer fragments S24, S26, S27, S28, S29, S31, and S32, respectively. Only the 
correspondence for fragment H agrees with that determined from the Horizons ephemeris, 
so we consider this one fragment to have a reasonably certain cross-identification with S28. 
The other suggestive identifications are indicated with '??' in Table [3) 

The lack of agreement between the ephemeris-predicted and Spitzer-observed positions 
is surprising, since the optical astrometry of small fragments that was used in the arc 
for the orbit solution stretches typically from 2006 Mar 24 through 2006 Apr 20, similar 
to the dates of the Spitzer observations. The estimated 3a error ellipse for the predicted 
positions varies depending on the length of the data arc; for example for fragments R, N, 
M, AB, BR, and Q, it is 28" x 0.2", 30" x 0.3", 64" x 3", 110" x 0.7", 330" x 0.8", and 
1517" x 1.5", respectively. Even taking into account these large uncertainties, the dispersion 
perpendicular to the orbit is significantly larger for the predicted (600") than observed (36") 
positions, and there should be more fragments with matches within 100" than observed. 

Why do the predicted positions not match the infrared fragments? Some possibilities 
include the following: (1) The ephemeris calculations are in error. This could arise due to 
the low-quality input data (inaccurate positions, which seem unlikely given the quality of 
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Fig. 8. — Zoom into the region around fragments 73P-M, N, and H. The top panel is an optical 
image (courtesy Peter Birtwhistle, Great Shefford Observatory West Berkshire, U.K.) taken 2006 
Apr 30 covering a 15' x 7' field of view with the identified fragments labeled in accordance with 
the astrometry reported to the Minor Planet Center. The bottom panel is a section of the Spitzer 
image taken on 2006 May 4 (during 12:49 to 15:50 UT), with a box showing the outline of the 
optical image of the top panel. Bold labels of the form S n with arrows correspond to the fragment 
list in Tab. [3j Circles with italicized labels indicate positions of fragments as predicted positions of 
fragments from ephemerides calculated by Horizons using orbit solutions based on optical, Earth- 
based astrometry. Vertical bars with labels at the bottom of the bottom panel are the approximate 
locations of fragments determined by P. Birtwhistle, based on their spacing in the optical image. 
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the input data and observers involved), or underestimated uncertainty when observations 
over a short arc from Earth are used to predict the locations of an object as viewed 
from a distance platform (given that Spitzer was 0.32 AU from Earth in May 2006). (2) 
Fragment identities were incorrectly matched to astrometry, e.g. different fragments were 
considered together as same fragment, or one fragment was artificially split into multiple 
short (uncertain) arcs. (3) Non-gravitational forces are significantly different from those of 
typical comets, with much higher perturbation perpendicular to the orbital plane (A 3 ) as 
well as within the plane. (4) The fragments do not last more than ~ 2 weeks, or are so 
highly variable in their outgassing that they are only briefly detectable, so that a different 
set of fragments was seen at each epoch. All of these explanations seem plausible for some 
fragments and may contribute to the lack of correspondence. We consider the last one in 
more detail here, since it leads to an important conclusion for the fragment longevity and 
evolution. 

The longevity of the fragments, and their brightness variations, can be surmised from 
amateur photometry compiled by the Minor Planet Center. The primary fragment C has 
a typical brightening and fading as a function of heliocentric and geocentric distances. 
Fragment B had a significant increase around 2006 Apr 1, followed by a decrease back 
to approximately the same trend that existed prior to the outburst. Fragment G had a 
dramatic increase in brightness around 2006 Apr 8, then faded rapidly — much more rapidly 
than a typical trend that would have peaked near close approach and stayed bright at 
perihelion, so that we suspect G has completely disrupted. Indeed, the HST images show 



G to be no more than a set of intermediate-sized fragments on 2006 Apr 19 (IWeaver et aL 



20061 ). The miscellaneous other fragments peaked in apparent brightness around close 
approach to Earth and then they too faded rapidly and were not seen after May 8. 

Looking back to the 2001 apparition, three components were reported: C and B are 
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the bright components from 1995 and again in 2006. But fragment E was for a short time 
brighter than B and then rapidly faded from view, decreasing in brightness after its sudden 
outburst while B and C continued their gradual increase toward perihelion. We suspect 
fragment E experienced a catastrophic disruption in October 2001 that made it briefly very 
bright (as an unresolved swarm of icy fragments) then rapidly fade from view. 

Thus it is plausible that 'lettered' fragments varied significantly in sublimation rate 
and sub-fragmentation, making their identifications problematic and their non-gravitational 
forces unpredictable. The orbits of the lettered fragments (other than C and B) are highly 
uncertain, and the infrared-detected fragment population (55 fragments observed together 
in one mosaic) differs significantly in distribution relative to the fragment C and B orbital 
planes. 
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Fig. 9. — Locations of fragments of 73P on 5 May 2006 as predicted using the JPL Horizons 
ephemerides (diamonds) and as detected with the Spitzer Space Telescope (asterisks). The 
horizontal axis indicates the distance behind fragment C along the line connecting fragments 
B and C, in units of the separation between B and C (Dbc)] the vertical axis indicates 
the distance of the fragment perpendicular to the line running through B and C, also in 
units of Dbc- The size of the diamonds depends on the magnitude from Horizons: those 
with no magnitude are smallest, those with magnitudes (H mag) greater than 20 are next, 
those with magnitudes less than 20 are next, and fragments B and C are shown with large 
symbols. The inset shows a zoom into a region particularly rich in fragments. The lack 
of correspondence between predicted and observed positions is evident, with the predicted 
fragments much more scattered perpendicular to the orbit than the observed fragments. The 
first astrometric observations included in the orbit calculations for all fragments was before 
the date of this Spitzer image. However, nearly all of the fragments were seen only briefly, 
and their latest reported astrometric observations included in the orbit determination were 
before the date of the Spitzer image. The fragments with astrometry reported after 2006 
May 4 are B, C, H, M, X, AI, AQ, BC, BN, BO, BP, and BR. A further 11 had astrometry 
up until 2006 May 1: G, L, N, R, AC, AS, BF, BI, BK, BQ, and BS. 
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Fig. 10. — Locations of fragments of 73P on 1 Apr 2006 as predicted using the JPL Hori- 
zons ephemerides (diamonds) and as detected with the Spitzer Space Telescope (asterisks). 
Only the positions of the very brightest fragments C, B, and G have corresponding infrared 
detections. The two fainter, infrared fragments, located along the orbit between B and C, 
do not match ephemeris-predicted positions of any fragments. Note that on 1 Apr, many 
of the fragments had not yet been optically reported (so they may not have existed or been 
activated). The fragments with astrometry reported before Apr 1 are B, C, H, J, K, L, N, 
M, V, U, P, Q, R, S, W, T, X, Y, and AU. Another 9 fragments were reported on 2 Apr, 
just a day after the first Spitzer image; these include Z, AA, AB, AC, AD, AE, AF, AG, and 
AV. The other 36 ephemeris positions plotted as diamonds are extrapolations back in time 
before the first reported position, so they may not be valid. 
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5. Properties of meteoroids 

The surface brightness of the trail on 2006 Apr 1 was measured from the IRAC 8 
fim and MIPS 24 /zm images. At 8 (24) /iin, the trail surface brightness in a region 
containing no significant fragments is 0.055 ± 0.007 (0.18 ± 0.015) MJy sr" 1 , the width on 
the sky is 148 (156)", and t he width correcting for dista nce is 6.9 (6.8) xlO 4 km. Applying 



extended emission (x 0.737, 



Reach and coauthors 



20051 ) and color (/1.06) corrections to 



the IRAC brightness; a color (/0.968) correction to the MIPS brightness; and scaling the 
IRAC brightness for the relative heliocentric distance difference to the MIPS observation 
(3 days later), the observed brightness ratio corresponds to a color temperature between 
the nominal wavelengths of 7.872 and 23.68 /im of T = 261 ±8 K. Scaling for convenience 
to 1 AU, assuming T oc r 5 yields T\ = 301 ±8 K. The trail is not detected at 70 /im in 
the 2006 Apr 1 image; the upper limit at 70 /im is 0.07 MJy sr -1 . Combining with the 
24 /im brightness, we can place a lower limit to the temperature of > 225 K, assuming 
greybody grains. Scaling for the heliocentric distance, we obtain a lower limit to the 
effective temperature at 1 AU of T\ > 259 K. This limit is consistent with the temperature 
inferred from the 8/im/24/im ratio. 

The debris trail is warmer than expected for dark, isothermal particles at 1 AU 
{Ti,iso = 278 K). The same result was 



60 (Ltm with IRAS by 



Sykes et al. 



bund for comet trails observed at 12, 25, and 
(120041 ) . who showed that the elevated temperature is 
consistent with that expected for rapidly-rotating, randomly-oriented spheres of low thermal 
conductivity (so they maintain a latitudinal temperature gradient). 

In principle, the elevated temperature could be due to physical optics, with the 
grains emitting less efficiently in the infrared than they absorb visible light; however, 
this possibility is definitively ruled out by the particle dynamics: particles small enough 
(< 20/im radius) for physical optics to be important feel such strong radiation pressure 
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that they cannot reside in the debris trail and instead are swept back and form the dust 
tail (with radiation pressure comparable to gravity). The distribution of meteoroids is 
diagnostic of their size and emission history. A detail ed analysis with comparison to 



2008|) . In 



theoretical models is presented in a companion paper (IVaubaillon and Reach! 
summary, the meteoroids in the debris trail are of order cm in size. The size distribution 
slope is similar for the meteoroids and in the size distribution around 100 /im to 1 mm such 
that larger meteoroids are much more abundant than would be guessed by taking the coma 
particle size distribution and extrapolating up to meteoroid sizes. 

The elevated infrared color temperature is more likely due to the particles being 
much larger than the thermal skin depth, such that a temperature gradient is maintained 
across the surface. For spinning bodies illuminated by sunlight, only those small enough 
for heat to be distributed from the day to night side more rapidly than the rotation 
period will be isothermal. For solid, rocky bodies the thermal skin depth is approximately 
L s = 1.7 x 10 3 P 5 fim where P is the spin period in seconds. A particle with diameter much 
smaller than the skin depth will be isothermal and can be treated as a dust grain, while a 
particle much larger than the skin depth should be treated like an asteroid. Particles with 
size comparable to the skin depth will have intermediate temperature distributions. To 



illustrate and bound the effects, we calculated temperatures o 



Dorschner et al. 



spher ical, isothermal grains 



19951 ) and compared them 



of amorphous olivine (optical properties from 
to the observed temperature and those expected for particles with temperature gradients 
Figure fTTl shows predicted temperature versus size. In the large-particle (asteroid) limit, 
the surface temperatures range from the subsolar temperature for a high-obliquity body 
(spin axis in the orbital plane; equivalent to no rotation or zero con ductivity) to the 



Harris et al. 



1998|) 



isothermal-latitude model applicable to a very fast-rotating body (cf. 
Since dynamical arguments rule out micron-sized grains, the models that can explain the 
observed temperature are those with low spin rates or moderate-to-low obliquity. A random 
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distribution of spin axes, as inferred by 



Sykes et al. 



(12004] ). is consistent with the observed 



temperature if the spin rate of mm-sized particles is less than 1 min. (Slower spins would 
yield higher t emper atures than observed. The spin states of meteoroids are not well known. 



Olsson-Steell (119871 ) calculated the paddle-wheel torque on particles from solar photons, 
finding a period of order 1 sec for mm-sized particles exposed to solar radiation for 1 month. 
The spin period is predicted to increase as the particle size squared, but it remains below 1 
min even for cm-sized particles. The predicted spin period for particles of mm to cm size, 
the size r ange inferred from the dyna mics modeling to contribute most of the debris trail 



emission (IVaubaillon and Reach 



20081 ). is consistent with the observed trail temperature. 
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0.1 1.0 10.0 100.0 1000.0 

radius (/xrn) 



Fig. 11. — Temperature versus size for spherical particles illuminated by the Sun at 1 AU. 
The solid curve was calculated using Mie theory for amorphous olivine. The other curves 
illustrate the behavior expected for particles with different surface temperature distributions. 
The asterisks are located at the skin depth appropriate for particle rotating once per second. 
The upper asterisk is for particles with obliquity (e) of 90°, i.e. spin axis in the orbital 
plane; the lower asterisk is for e = 0, i.e. spin axis perpendicular to the orbital plane. 
The diamonds are located at the skin depth appropriate for particles rotating once per 
minute. Dashed and dotted lines are simple interpolations between the isothermal limit (for 
amorphous olivine) and the temperatures of particles with skin depths for 1 sec and 1 min 
periods, respectively. The 2a confidence interval for the observed temperature is shaded as 
a grey rectangle, truncated at 100 /zm particle size based on dynamical constraints. 
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6. Dynamics and survival of fragments and meteoroids 

We now address the dynamics of the fragments. We differentiate between the primary 
(B and C) fragments, the large fragments (lettered by the Minor Planet Center and visible 
as a distinct miniature comet in the Spitzer images), the small fragments (visible only 
as individual entities in the Hubble images and otherwise as the lumpy inner coma of 
ground-based optical images), and the meteoroids or debris stream (visible as extended 
emission along the orbit only in the Spitzer image). 



6.1. Sublimation and survival 

The fragments are undoubtedly composed of ice and refractory material. A fraction of 
the surface /j ce is exposed, and the sublimation rate of the ice is defined 

pT x dM _ 2 
Z = — — — — mol cm s (2) 

dAdt v ; 

where dM is the mass of ice produced per unit surface area dA and time dt, and \i is the 
molecular weight of the ice (presumably H 2 0). The time that a fragment can survive before 
its ice is sublimated is 

t. = — = (3) 
fiZA 3/j,Zf ice 

which in physical units becomes 



*. = ( ( A) (1*2£) Z-i days, (4) 

\gcm 6 J \cm/ \ n J 

where Z\$ is the sublimation rate in 10 16 cm -2 s -2 , R is the fragment radius, and p is 
the fragment mass density. The sublimation rate depends on the particle temperature. 



A surface o 



(IDelsemme 



pure water ice illuminated at 1 AU from the Sun will have Z 16 = 1.7 x 10 2 



2000). The sublimation rate at the surface of a comet depends on its 



temperature distribution, spin rate, conductivity, and composition, with estimated values of 
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order Z 16 ~ 10-30 fjWhipple and Huenberlll976l ). Small ice particles are colder than a large 



slab and hence su 



jlimate more slowly, with Z 1fi ~ 0.2 for R = 0.1 cm and Z 16 ~ 7 x 10 4 



1975 



for R = 0.01 cm (IPatashnick and Rupprecht 
be warmer, with Z 16 ~ 10 3 for R = 0.5 cm dMukai 



. Het erogeneous particles may however 



19861 ). 



The small fragments must survive at least over the 2 HST orbits during which they 
were observed to move rapidly but could be traced as distinct individuals; setting t s > 2 
orbits requires R > 400/j ce cm. This condition is likely to apply since the inferred size 
of the small HST fragments is much larger; however, the HST magnitudes sizes could be 
overestimated if they are dominated by comae. The HST fragments that disappeared 
between observations (or changed into diffuse clouds) probably disaggregated into debris, 
visible in the images due to dust scattering sunlight, rather than sublimating into vapor. 

The debris trail particles must survive years in order to occupy their observed 
locations along the projected orbit; this requires that the debris trail particles have radii 
R > O.QficeZiQcm. In terms of lifetime, small amounts of ice could be present; however, the 
dynamics arguments below show that fi ce must be small in the trail particles. 

The basic model for cometary grain production has the expanding gas from sublimating 
ice near the surface carry off solid materials that are entrained if the drag force is greater 
than the gravitational pull of the nucleus: 



1 2 GMm 

-pvZ.a > 



(5) 



where v is the expansion velocity of the vapors from the sublimating ice, a and m are 
the particle radius and mass, Rn and M are the nucleus radius and mass, and G is the 
gravitational force constant. The maximum size that can be levitated is then 

9fivZ 



(6) 



32itGpnpRn 

where a and p are the grain radius and density, and p^ is the nucleus density. In physical 



-46 - 



units, for a nucleus with sublimation rate Ziq = 20, 

dmax ~ 10 (^^j Cm > ( 7 ) 

for p = pN = 1 g cm" 3 . Since 73P has perihelion g = 0.94 AU and the nuclear ra dius 



in 1994, before the recent breakup events, is Rn < 1.1 km (IBoehnhardt et al 



19991 ). the 



largest particles that can be levitated by gas drag are of order 10 cm in size. 

After decoupling from the gas, the particle's trajectory is determined by the velocity 
that has been imparted to it by gas drag, v e j, and the ratio of radiation pressure force to 
solar gravity, 

P = Frad/Fgrav = 0.57Q pr p _1 (1 /an/a) (8) 
for spheric al particles, where Q pr is the a coupling efficiency between the particle and solar 



radiation (IBurns et al. 



19791 ); for particles much larger than the wavelength of visible light, 
i.e. o> 1 /jiii, we can set Q pr = 1. For fragmentation or disaggregation of a nucleus (or 
meteoroid), the treatment is analogous to that of debris, with the velocity imparted to the 
fragment (due to whatever process caused the fragmentation) taken to be v e j. 

For small objects, and additional term must be added to account for the 'rocket' effect 
due to expansion of vapors from sublimation of ice on its surface. To explore this effect, we 
define 

^ Frocket 3flZVi ce fi ce 

Fgrav 4GM G pa 

where v i ce is the velocity of expansion of the sublimating ice. In physical units, 

Using this parameterization, the rocket effect a is analogous to radiation pressure j3: both 
represent forces that point away from the Sun and are fractions of solar gravity. It is 
furthermore useful to make the ejection velocity dimensionless through v = v e j/vorb, where 
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v or b is the orbital speed at the time of particle production (normally, close to perihelion). 
The parameters (3, a, v allow a dimensionless means to determine the dominant celestial 
mechanical perturbation for a particle. (Note that unlike the radiation pressure parameter, 
(3, the rocket parameter, a, depends on heliocentric distance. This is due to the fact that 
while insolation and radiation pressure decrease with distance from the Sun as r~ 2 , the 
sublimation rate has a steeper dependence taking into account the need for the ice to be 
heated above its latent heat before it will sublimate. Furthermore note that we assume 
the outgassing is in the direction of the Sun, so that this parameterization neglects the the 
rocket force due to the component of the outgassing that is nonradial.) 
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^~*^© 10 s km 




Fig. 12. — Illustration of dynamics for (a) dry particles and (b) icy particles. The arrows 
at the top indicate the direction of the velocity vector (v) and the antisolar direction (— 0) 
projected onto the observed sky plane. The dynamics of dry particles is determined only by 
the ratio of radiation pressure to solar gravity, (3. Particles in the anti-solar direction, such 
as observed in the tails of the 73P-B, 73P-C, etc fragments, have (3 > 10 _1 corresponding 
to sizes < 10 fim. Particles along the debris trail (which closely follows the projected orbit) 
must have (3 < 10 -4 , corresponding to sizes > 1 cm. The dynamics of icy particles is 
determined by the ratio of rocket force (from sublimating ices) to solar gravity, a, as well 
as (3. The particle illustrated in panel (b) has a radius 1 mm, (3 = 10~ 3 and begins with 
fice = 0.0017 so a = 10 _1 . Since initially a ^> (3 > u, the particle follows the syndne for 
a = 10 _1 . The particle becomes devolatilized (no more ice exposed to sunlight) at time r^, so 
a — > 0, and the remaining, effectively dry particle then follows a (3 — 10~ 3 trajectory. Debris 
trail particles follow the projected orbit closely and therefore must have — > 0. The small 
fragments in the HST images follow the tail-ward (large a) direction, then either disintegrate 
or follow low-/? trajectories so far from the parent nucleus that they are effectively lost from 
the system. 
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6.2. Dynamical classes of fragments 

The trajectories determined by the combination of v and a + (3 can be divided into 
classes. Large fragments have a + (3 ~ 0; they remain close to the projected orbit, with 
trajectories determined only by v. (Gravitational perturbations by the planets are of 
course important on long timescales. Here we consider only the 3 orbits during which 
the fragmentation events transpired and during which no close encounters with planets 
occurred.) The separati on velocities fo r fragments of split comets have been deduced in 



the range of 1-2 m s _1 ( ISekaninal Il978l ) . corresponding to v ~ 10 -5 . For comparison, the 
radiation pressure has (3 < 10 -5 for particles with radii a > 6p _1 cm. Intermediate- sized, 
icy fragments appear in clusters around large fragments (perhaps being their principal 
building blocks) and have dynamics, relative to their local center of mass, dominated by 
a, the rocket effect; the size scale is of order 10 m. These fragments are discussed in more 
detail in a paper describing the HST results; the dynamical modeling and discussion are 
included in this paper in order to contrast these fragments from the others. Debris trail 
particles are smaller and lie along the orbit of the parent comet; they have evidently a <C (3 
so that their dynamics is dominated by radiation pressure. In the following subsections we 
examine the dynamics of each type of fragment in more detail. 



6.3. Separation of large fragments 

The small perturbation caused by the radial and in-orbit components, v», of the 
ejection velocity causes large fragments with a + (3 ~ to spread along the orbit of the 
nucleus, while the perpendicular component, v±, causes fragments to disperse perpendicular 
to the orbit. We can therefore use the observed locations of the fragments perpendicular 
to the orbit to estimate vn and v±. A natural but wrong way to compute these quantities 
is to simply divide the distance between the fragment and the comet by the time since its 
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ejection; however, care must be taken here. A fragment ejected with the velocity v e j at a 
point P of the cometary orbit will have a slightly different inclination % and semi-major 
axis a than the main body. Considering Keplerian motion, this implies that the fragment 
will always pass by point P at each subsequent passage. As a consequence, the physical 
separation between the fragment and the main body will vary with the true anomaly. We 
assume all fragments were produced in the 1995 breakup event, and we determine vn and 
v± by matching the location of the fragments by iteration. Because v\\ affects a, it has a 
strong influence on the location of the fragments within the orbit plane. Nearly all the 
large fragments are located closer to the cometary orbital plane, and their inferred v» ~ 1 
m s -1 . Interestingly, their inferred v± is similar to v\\, suggesting the ejection from the 
parent comet was roughly isotropic. 

One fragment, illustrated in the inset at the top of Fig. 3b, is located far from the 
comet orbital plane. Its distance from fragment B in the plane is 1.9° (or 1.8 x 10 6 km), 
and it is 4.1' (or 6.6 x 10 4 km) away from the orbit plane. Assuming that this fragment 
was ejected during the 1995 breakup event, we find vn ~ 1.4 m s" 1 and v± ~ 16 m s _1 . The 
value of v± is very high, relative to v\\ for this fragment and inferred ejection velocities of 
other comet fragments. (All the other fragments are located much closer to the cometary 
orbital plane, with corresponding v± is less than 2 m s _1 .) 

The location of a fragment can be influenced by non-gravitational force from 
asymmetric outgassing (a ^ 0). Translated into the Marsden (1961) non-gravitational 
force parameter, we find A 3 <, 5.0 x 10~ 9 for all fragments except the one discussed in the 
previous paragraph. This value is higher than typical Jupiter-family comets, where the 
median absolute value is \A 3 \ = 1 x 10~ 9 . The fragment located furthest from the orbital 
plane has Ai = 8.9 x 10~ 9 and A% = —2.0 10~ 7 . The later value is one order of magnitude 
higher than the highest value recorded for cometary nuclei, namely comet C/1996 Nl 
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Brewington (A 3 = -1.3 1(T 8 ; P. Rocher, IMCCE, cometary note 020s[j). however given 
the size of the fragment one expects that the outgassing process will have a strong influence 
on its dynamics as discussed above for the mini-fragments in the HST images. 

The reality is that a combination of ejection velocity and non-gravitational forces 
determine the fragment locations. Our findings show that the fragments observed by Spitzer 
have high but reasonable ejection speeds. One of the fragments must have been released at 
very high speed perpendicular to the orbit. It must be stressed that none of these fragments 
could have been ejected in 2006 — otherwise, the ejection velocity or Ai values would be far 
higher than the values reported here. 



6.4. Rocket effect and HST fragments 

The rocket effect is more important than initial ejection velocity (a > v) for particles 
with ap < 3.Qfi ce v e j km. That non-gravitational rocket forces are significant for objects 
the size of cometary nuclei is well-known, but the ice-sorting effect on fragments is a new 
application. The effect is analogous to radiation pressure in direction, but is approximately 
G00Zi6fi ce times larger, so the rocket effect and radiation pressure are equivalent for 
fi ce = 10 -2 . The small fragments seen in the HST images are estimated to have sizes 
of order 5 m. Some of them are visible in 3 consecutive HST orbits, so they survive at 
least at least 200 min. Their apparent motion on the s ky is ~ 100 m s" 1 i n the antisolar 



direction — the same direction as the dust and gas tails (IWeaver et al\ 120061 ) . The direction 
of motion clearly indicates a + (3 > 10 -2 , and the size indicates /?Ca. Combining the size 
and a, we solve for the fraction of the surface of the small fragments that must be covered 



2 http: / /www. imcce.fr/en/ephemerides/donnees/comets/FICH/CIA0209.php 
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by ice 

/ a \ /km s _1 . 

/ice > 0.8 (— ) , (11) 



s5 Hi/ V ^ice 

so the small fragments must be very icy in order to move in the tailward direction. It is 
possible that the sizes of the fragments are overestimated from the HST images due to their 
bright comae, but the presence of those comae demonstrate that the fragments are very 
active and icy; thus fi ce is unlikely to be smaller than 0.1. 

A significant difference between the rocket and radiation pressure effects is that the 
former may be time-dependent, if the fraction of exposed ice changes with time. We discuss 
this in more detail in the next subsection with regard to the meteoroid-sized particles. The 
fate of the icy fragments is to either (1) continue along high-a trajectories, (2) become 
'extinct' and transfer to orbits roughly similar to the parent comet, or (3) disaggregate 
into clouds of dust and gas. Even in cases (1) and (2), it is clear that the small fragments 
are lost from the cometary system, in the sense that their orbits are so different (and the 
objects are so small) that they would be unlikely to be detected again and even if so their 
orbits have already been sufficiently perturbed from that of the parent comet that they will 
be many degrees from that comet upon return and would not be noticed in observations 
centered on the comet. Thus a general result is that fragments in the ~ 5 m range are 
a form of missing mass that could only have been detected thanks to the high angular 
resolution of HST and the close approach of 73P to the Earth. Observations of another 
fragmenting comet LINEAR indicate d that the detectab le fragments and dust do not add 



up to the comet's pre-breakup mass (jWeaver et al. 



20011 ) . It now appears that a significant 



portion of the cometary fragmentation may be in the form of icy fragments in the 5-m size 
range. 

Another aspect of the small fragments seen with HST is that some of them disappear 
in between successive images, and their are some fuzzy patches without nuclei that move 
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away from the comet. This fragmentation and disaggregation appears to be analogous to 
that of the main fragments, except that for the smaller ones the there remains only an a 
cloud of vapor. The material strength of the fragments is therefore probably small. Based 
on the tidal splitting of Shoemaker-Levy 9, the strength of the nucleus was estimated to 



be ~ 10 dyne cm 2 (lAsphaug and Benz 



1996). If the small fragments have strengths of 



this order, then they are even more likely to fragment and disintegrate because they have 
far less gravity than the nucleus and cannot support themselves as 'rubble piles' once they 
have been damaged by collision or thermal stress. 



6.5. Radiation pressure and debris trail 



The debris trail in Figures [2] and |4] clearly shows the presence of long-lived debris 
within th e comet's orbit. We can con strain a + (3 < 10~ 5 from the location and width of 



the trail (jVaubaillon and Reach 



2008) . If the part icles were levitated from the nucleus, 



3), 



20071 ). then they have radii < 10 cm and 



as in normal cometary activity (IReach et al. 
(3 > 6 x 10~ 6 . Thus only the largest particles that can be levitated from the nucleus reside 
within the debris trail, and there is no room for these particles to have significant rocket 
effect. 

Larger particles could be present in the debris trail if they were released from the 
nucleus by another mechanism like nuclear fragmentation. They must still achieve escape 
velocity from the nucleus, v esc ~ l(/2jv/km) m s , which corresponds to a velocity 
perturbation parameter v > 2 x 10~ 7 . Meteoroids of order 1 m size could be present in the 
debris trail with small amounts of ice as long as a < 10 -5 . Referring to Eq. [101 we can then 
constrain the sublimation rate Z w f ice < 0.03. The size range discussed in this paragraph 
defines a transition between the icy fragments that follow the dust tail, like those seen 
in the HST images of 73P, and smaller and lower-sublimation-rate bodies. Note that the 
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observed surface brightness of debris trails is unlikely to be produced by m-sized chunks; 
the required mass of debris would be too high. Furthermore, meteor showers are due to 
~mm-sized meteoroids. The m-sized fragments would be rare denizens of the trails. 

As volatiles disappear from the surface of a debris trail particle, its interior may 
eventually become either dry or shielded by refractory material remaining at the edge 
of the particle. The situation is somewhat analogous to that of the comet as a whole, 
where a 'rubble mantle' of particles that cannot be levitated eventually shields the interior 
from sunlight and reduces the active area so fi ce <C 1 and the activity and its resulting 
non-gravitational force terminate. Cometary debris particles have negligible gravity, but 
their material strength may be significant. 

If the refractory part of a particle had no material strength at all, e.g. if it were an 
agglomeration of very small particles held together only by the ice, then the particles would 
disaggregate while they are sublimating. The tiny subunits would flow away at the same 
time as the gas, and the particle would simply shrink until only one subunit is left. If one 
believes this model for the cometary debris constitution, then a is always large, it dominates 



radiation pressure, and the lifetime is set by Equation [H Cometary debris trail partic 



survive for > 10 yr (two orbital revolutions ISykes and Walker 



1992 



Reach et al. 



20071 ). 



cs 



Averaging over a typical short-period comet orbit with 5 yr period, {Z\$) ~ 0.2Zi6,i, where 
Zi6,i is the sublimation rate at 1 AU. For particles with Zi§^ ~ 1 to survive for > 10 yr, 
they must be larger than 14 cm in size. Such particles could not have been levitated from 
the nucleus by gas drag, and the mass of such particles required to explain the observed, 
smooth and bright, debris trails would be orders of magnitude greater than that of the 
nucleus from which they origi nate. Meteoroids producing known meteor showers survive 



for > 10 3 yr (jjenniskens 



2006); only bodies much larger than 1 m could survive continuous 



sublimation so long. 
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That the dynamics of the debris trails and the small fragments are distinctly different 
indicates that their internal structures are different. The debris trail particles must be 
dry. Based only on the location of the particles, the size a < 10 cm, and a < 10~ 5 , so 
fi ce < 2 x 10~ 5 . The survival of the debris trail particles requires that they either originated 
from a completely dry portion of the nucleus or they have become devolatilized almost 
immediately after separation from the nucleus. If the particles originated from a completely 
dry portion of the comet, then we must explain how they were ejected from the comet. 
Possibly the particles were from the mantle, and they were lifted from the surface upon 
activation of a patch of buried ice. Alternatively, the particles could be from a part of the 
nucleus adjacent to an active region, and they fell into or were otherwise entrained in the 
outflow. Both of these scenarios and the subsequent dynamics of the particles are consistent 
with the interplanetary dust particles of cometary origin having anhydrous mineralogy 



([Messenger et al 



empel 1 and 81 P /Wild 2 appearin 



2005 



Brownlee 



20041 ). If 



20061 ) and the surfaces of 9P/ r 
predominantly 'dry' in recent spacecraft fly-bys flA'Hearn et al. 
the particles did maintain some volatile compounds, then in the process of devolatilization 
upon release from the surface and exposure to the full ultraviolet radiation of the Sun, the 
refractory parts of the particle must have retained strength. That is, they could not be 
composed of very fine (< 10 /im) grains that are held together only by volatiles acting as 
a sort of 'glue.' Thus the debris trail particles (and meteoroids of the size that produce 
meteors) must have some material strength. 

This inference, based only on dynamics, is compatible with knowledge of the properties 
of meteors and interplanetary dust particles. Of some importance for understanding the 
nature of comets and interpreting collected, or remotely observed, small (< 100 /im radius) 
interplanetary particles in comparison to the larger (> 1 mm) particles that dominate 
the total mass loss, is whether the larger particles are amalgams of the smaller particles, 
whether the small particles are amalgams of fine grains, and how those subunits are bound 
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together. Meteors are known to fragment in the Earth's upper atmosphere, which suggests 
they are not very strong, but in fact the impact onto the Earth's atmosphere is an energetic 
event and still allows for significant streng th, estimated in the range 10 4 dyn e cm 2 with 



the weakest ones at ~ 4 x 10 2 dyne cm 2 (jTrigo-Rodriguez and Llorca 



2006J). The stress 



on the particle's refractory structure can be roughly estimated from the sublimation rate 
ZvfirriH ~ 5 dyne cm -2 for sublimation at 1 AU from the Sun. The sublimation rate will 
scale as r -2 , so particles with the strength of weak meteors would explode if they were 
released within 0.1 AU of the Sun. (The existence of Geminid meteoroids, with perihelion 
0.13 AU, argue for a high strength for bodies at ~mm sizes.) But for a short-period comet 
like 73P (or even the very-low-perihelion short-period comet 2P, with q = 0.339 AU), such 
close approaches to the Sun to not apply and particles have sufficient strength to survive 
sublimation of their embedded ices. 

Meteor showers have been associated with comets through similarity of their derived 
orbits. The meteoroid orbits are usually different from the present orbit of the comet, 
and are instead associated with past orbits. The primary acting agent in the celestial 
mechanics of meteoroid streams is gravitational perturbations by Jupiter. Even in absence 
of a close encounter with Jupiter, the orbit of the comet slightly changes at each return. 
The meteoroids are ejected on orbits similar but independent from the comet. Because of 
the ejection velocity and the radiation pressure, their orbital semimajor axis is typically 
greater than that of the comet. As a consequence the Jovian perturbations are not the 
same at each point of the stream, and the result is the deformation of the whole structure. 
As long as the perturbation are small enough it is still possible to associate a trail with its 
parent body. However, the knowledge of the passage of ejection requires a full calculation 
of the evolution of the stream. The success of s uch calculations in pre dicting the peaks 



of the Leonid showers confirm the basic theory (jVaubaillon et al 



20051 ). If the meteors 



are ejected from the comet with some surface ice with f ice > 10 2 , then their initial orbits 
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will be affected. Simple treatments of the physics of separation from the nucleus consider 
one parameter, the ejection velocity, which includes effects of gas drag, decoupling, and by 
extension effects such as sublimation of retained ice. The total work done on the particle is 
the rocket force integrated over the path until devolatilization, which we can characterize 
with the ratio 

w = — — = ar dv /P (12) 

" grav-t 

where P is the orbital period and r dv is the lifetime of ice on the surface. Constraints 
placed traditionally upon the ejection velocity may be alternatively interpreted as a limit 
on w with order-of-magnitude w ~ v e j/v or b = v. A limit of v e j < 10 m s _1 corresponds to a 
limit of w < 0.2 for typical orbital velocities near perihelion. For an orbital period of 5 yr, 
we then find a limit ar^ < 1 yr. For cm-sized particles, we then find a limit r^/j ce < 0.2 
yr. Thus if meteoroids are very icy (fi ce > 0.1), they must become devolatilized within 2 yr. 
Conversely, the ice that can remain in meteoroids that are 10 2 yr old is fi ce < 2 x 10 -3 . 



7. Gas in the coma 

The IRAC images at 4.5 /im have a morphology distinctly different from that at the 
other wavelengths. The spectral energy distribution (Fig. IT31 clearly shows that there are 
locations in the coma of 73P-B with a distinct excess of 4.5 /im emission with respect to 
the continuum that dominates most infrared wavelengths. Figure [13] shows the spectral 
energy distributions of a sunward portion of the coma and and anti-sunward portion of 
the dust tail. The continuum brightness due to dust is a combination of thermal emission 
and scattered light. The scattered light has negligible contribution except possibly at the 
shortest wavelength: if all the 3.6 /im brightness were scattered, the estimated scattered 
light in the other bands (i.e. the 3.6 /im brightness scaled by the solar spectrum) is <10% 
To fit the thermal emission for the dust tail, we use the 5.8 and 8 /im brightness, after 
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color correc tions are made for the wid e IRAC bandpasses as appropriate for a 300 K 



blackbody (IReach and coauthors 



20051 ). to measure a color temperature of 288 K for the 
tail. The scattered light brightness required to explain the 3.6 /im emission above the 288 
K blackbody emission amounts to 50% of the 3.6 jttm brightness and 6% of the 4.5 /im 
brightness. 

For the forward coma, we remove the combined thermal emission and scattered light 
assuming the spectrum of the dust is the same as in the tail. While the dust temperature 
is likely to vary somewhat within the coma due to particle size changes, the variations are 
evidently not large, because the tail spectral energy distribution falls very close to the coma 
brightness at 5.8 and 8 /im. For the region used in Figure [TBI the excess 4.5 /im brightness 
amounts to 0.75 MJy sr _1 within the IRAC band. We attribute the excess brightness to 
line emission. Based on the spectrum of 103P/Hartley 2 (a Jupiter-family comet similar to 



jandpasses, by 



ar, is the 



1999|). The 



73P) obtained with ISO, the brightest feature in all the IRAC 
CO2 V3 band at 4.26 /im, with a weaker contribution from CO (iColangeli et al\ 
brightness of the excess emission attributed to gas features in the forward coma of 73P-B, 
at 5000 km from the nucleus, is 1.1 x 10~ 4 erg s _1 cm -2 sr _1 . 

The distribution of line emission in the 4.5 /zm band (after removing thermal emission 
using 288 K blackbody scaled to match the 5.8 and 8 /im bands and scattered sunlight 
using a solar spectrum scaled to match 50% of the 3.6 /im band) is shown in Figure [TH It 
is already evident from the 4.5 /im image alone that there is an extra spatial component 
located in the sunward direction. The subtracted image shows the distribution of this excess 
is a very wide elliptical fan-like region primarily forward of the nucleus. The full-width at 
half maximum brightness is 36" (17,000 km) and the emission is detected out to 79" (37,000 
km) in the sunward direction. 

The radial profile of the CO2 roughly follows a shape, where p is the separation 
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wavelength (/xm) 



Fig. 13. — Spectral energy distribution from 3.6-8 /jm obtained from IRAC photometry of 
two locations near 73P-B. The upper curve is the forward coma brightness, averaged over 
a circular region with radius 7" centered 11" (5000 km) in the sunward direction from the 
nucleus; it is scaled by a factor of 5 for clarity in this plot. The lower curve is the tail 
brightness, averaged over a circular region with radius 30" centered 197" (92,000 km) in the 
anti-sunward direction from the nucleus. The solid line simply connects the data points. The 
dashed line is a blackbody fit (T = 288 K) to the tail brightness, shown both through the 
tail brightness measurements and again through the forward coma brightness (after scaling 
to match the 5.8 and 8 fxm points). The dotted line is the tail brightness, scaled upward to 
match the 5.8 and 8 fim forward coma brightness. There is a clear excess of emission in the 
4.5 /im band in the forward coma, with respect to the dust spectrum of the tail. 



Fig. 14. — IRAC images of the coma and part of the dust tail of fragment B, 12' x 7' in 
angular size: top left: 3.6 /im band, top right: 4.5 /im band, bottom left: 5.8 /im band, bottom 
right: excess 4.5 /im band after subtracting scaled 3.6 and 5.8 /im bands. The 3.6 /im and 
5.8 fim images are dominated by dust scattering and thermal emission, respectively. The 
excess 4.5 /im is due to CO and C0 2 gas emission. 
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from the nucleus on the plane of the sky, as expected if the CO2 molecules have a r 2 spatial 
density profile and observations are within the disso ciation distance . The dissociation 



1994 ). so the dissociation 



timescale for CO2 is 5.8 days at 1 AU from the Sun (jCrovisierl 
scale length corresponds to 6' from the nucleus in 2006 May; that the dissociation length is 
further than the distance out to which we measured the radial profile indicates dissociation 
is not expected to dominate. A profile would apply for outflowing molecules in steady 
state produced from the nucleus, as opposed to an extended source, e.g. sublimation from 
surface of grains or production by photo-dissociation of larger molecules. That the observed 
radial profile is consistent with a nuclear origin is consistent with CO2 dominating the line 



emission in the 4.5 /im band, consid ering that CO has been 



source in many cases, as reviewed by 



Bockelee-Morvan et al. 



bund to have an extended 



(120041 ) . Integrating the profile 



out to 48" radius leads to a flux 8.6 x 1 12 erg cm 2 s 1 . T he flux is a factor of 3 lower 



than that of the very active comet 103P (jColangeli et al. 



1999). so using the same excitation 



scaling, the production rate Q(C02) ~ 1 x 10 27 s _1 . This flux must be treated with caution 
since the dust continuum subtraction is inaccurate near the nucleus where the emission is 
brightest; the total flux also scales linearly with the radius out to which one integrates. 

For fragment C, the spectral energy distribution of the tail and coma are similar to 
that of fragment B. There is also a 4.5 fj,m excess from the sunward side of the coma, with 
similar morphology, though it is lower in amplitude than that for fragment B. Applying the 
same spectral and spatial fit as described above, we find a flux ~ 3.0 x 10 -12 erg cm" 2 s -1 , 
and an estimated production rate Q{C02) ~ 7 x 10 26 s -1 . 



Kobayashi et al. 



(120071) on 



The H2O production rate of fragment B was measured by 
2006 May 7, when the comet was at 1.03 AU from the Sun, to be 3.5 x 10 28 s -1 . If both the 
CO2 and H2O production rates scale as r~ 2 , then the ratio of production rat es at 1 AU is 
Qco 2 /Qh 2 o ~ 0.05. The H 2 production rate of fragment C was measured by lDiSanti et al. 
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( 120071 ) on 2006 May 28, when the comet was at 0.95 AU from the Sun, to be 1.3 x 10 28 
s _1 . If both the CO2 and H 2 production rates scale as r -2 , then the ratio of production 
rates at 1 AU is Qco 2 /Qh 2 o ~ 0.11. Evidently C0 2 is highly abundant in the 73P-B and C 
ices. The abundant CO2 in 73P could be due to its recent exposure to sunlight since the 
splitting event in 1995 (and subsequent splittings for fragment B), or it could be the typical 
cometary composition. For 103P/Hartley 2, which has not split re cently and was observ ed 



19991) 



at 0.84 AU from the Sun, the C02/H 2 abundance ratio is 0.08 (jColangeli et al. 
comparable to the 73P fragments. 

For comparison, the CO producti on rate was found to be exceptionally low (relative to 



H 2 0) for 73P-C: Qco/Qh 2 o = 0.005 (IDiSanti et al\ 120071 ) . Combining the near-infrared 



spectral data with our 4.5 /jm result yields an abundance ratio Qco 2 /Qco ~ 20 making 
C0 2 the second-most-abundant (by f ar) molecule in the ice exposed on 73P-C. For comet 



Hale-Bopp at 2.93 AU from the Sun, 



Crovisier et al. 



( 119971 ) measured production rates 



with relative proportions for H_20:CO:C02 of 100:22:70. The production rate of CO2 was 
very high for the Hale-Bopp observation (even higher than we find for 73P), but that is 
mostly due to the large heliocentric distance for the Hale-Bopp spectra — the H2O remains 
solid and does not sub l imate rapidly at 2.9 AU. After extrapolating inward to 1 AU, 



Bockelee-Morvan et al. 



(120041 ) estimated production rates for Hale-Bopp with proportions 
H 2 0:CO:C02 of 100:23:6. Thus the CO2 production rate relative to H 2 for 73P-C appears 
to be within the range of other comets, and the CO is anomalously low in the 73P-C ices. 
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8. Comparison to other comets 
8.1. Icy fragments 



For C/1996 B2 Hyakutake, iDesvoivres et al\ (120001 ) identified and measured 



trajectories of 7 distinct condensations in the coma. The separations of the fragments 
increase quadratically with time, clearly indicating acceleration due to non-gravitational 
forces and allowing an accurate measure of the acceleration and the date of splitting from 
the main nucleus. Taking the observed separation and time since splitting, t, we can 
determine the a due to the rocket effect as defined above, 

a = ¥gm~ & ' (13) 

where d is the separation from the nucleus corrected for the projection of the solar direction 
on the sky. The 7 observed condensations have 1.2 x 10~ 2 > a > 5.4 x 10~ 4 . These values 
are in the same range as those required to match the fragments observed with HST and 



Subaru for 73P. 



Desvoivres et al 



( 120001 ) present a physical model for the non-gravitational 
force in terms of ice sublimation and derive from it a range of fragment sizes (times the 
density, p in g cm -3 ), of 1.4 < pa < 24 m. Again, these values are comparable to what we 
derived from the dynamics and the photometry of the 73P-B fragments in the HST images. 

Future high angular resolution imaging of close-approaching comets could reveal 
how common is ice fragment production and how important a mass loss mechanism 
is in general. It is evident that production of large icy fragments is not exclusively a 
property of long or short-period comets, given their detection from C/Hyakutake and 
73P/Schwassmann-Wachmann 3. It is clear that production of icy fragments is a significant 
mass loss mechanism, when it is happening. For Hyakutake, if we sum only the mass of 
the 7 detected fragments, and assume these are the only ones produced over a month of 
observation, then the mass production rate is 1.2 x 10 5 g s -1 . Of course, it is unlikely that 
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Hyakutake is fragmenting at the same rate throughout every orbit; if the 7 fragments from 
1996 were the only fragments produced over an entire orbit (10 5 yr), then the production 
rate is a mere 0.1 g s _1 . 

For 73P-B, from the HST or Subaru images, we estimate at least 20 fragments with 
sizes of order 10 m produced over a period less than one week, yielding 1.4 x 10 5 g s _1 mass 
production rate. During the Giot to flyby of comet 26P /Grigg-Skjellerup in 1992, a fragment 



of order 10-100 m was detected (IMcBride et al. 



19971 ). meaning three very different types 



of comets have intermediate-sized fragments that are too large to have been levitated from 
the surface by gas drag, and, even if they survive long enough, too small and too far from 
the nucleus to ever be found again as separate comets. It is not known over what fraction 
of the orbit such fragmentation occurs, but it does not seem to be rare. 



8.2. Debris trail and swarm 



Debri s trails have been detect ed for more t han 30 comets usi ng infrared observations 



by IRAS flSvkes and Walker 



19921 ) and Spitzer (IReach et al. 



20071 ). A debris trail due to 



73P is clearly evident in 



IRAS flSvkes and Walker 



Spitz er images. T he debris trail fro m this comet was evident in 



1992|) and COBE flLisse et al 



19981 ) data before the 1995 splitting 



event. Only infrared images have revealed the debris trail; it has not yet been detected 
optically despite extensive observations with a wide range of telescopes in 2006. Debris 
trails are due to particles in the mm to cm size range, and the orbit-averaged mass loss rate 
is of order 10 4 g s -1 , comparable to the mass production rate of H2O. Averaged over entire 
orbits, the mass loss rate in the form of debris appears to be much larger than the mass loss 
rate in the form of fragments. 



The properties of the diffuse infrared emission in the Spitzer images indicate that the 
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particles are somewhat different in origin than typical debris trails. While the surface 
brightness along the orbit and ahead of fragment C, in the 2006 Apr, 2006 May, and 
2007 Jan images, is typical of debris trails, the brightness following fragment C, and in 
particular the significant enhancement of diffuse emission behind fragment B in 2007 Jan, 
is highly unusual and cannot be explained by the same origin as typical cometary debris. 
For a quiescent comet, the production rate increases towar d perihelion, modu lo active 
region seasonal effects such as are important for 2P/Encke (IReach et a/.ll2000l ). The mm 
and cm-sized debris gradually separate from the comet; they remain close to the nucleus 
during the comet's perihelion passage (due to the small ejection velocity and small effect 
of radiat ion pressure). T 



streams (IJenniskens 



'hese large particles form debris trails, or equivalently, meteoriod 
20061 ). that follow the comet's orbit on the subsequent revolutions. In 
this model, the brightness profile along the trail is relatively smooth, gradually decreasing 
behind the nucleus and occupying more and more of the comet's orbit until a gravitational 
perturbation (normally by Jupiter) changes the comet's orbit and forces the comet to start 
anew in building its meteoroid stream. 

For comet 73P, the brightness profile through the 2007 Jan image is neither monotonic 
nor smooth. Instead, as shown in Figure [15j there is an elevated brightness behind fragment 
C that remains fairly constant, then a significant rise behind fragment B followed by a 
gradual decline. The broad 'bump' in the surface brightness profile behind fragment B 
is centered around -1.4 times the separation between fragments B and C (Dbc), relative 
to fragment C; the width of the bump is approximately 0.3 (FWHM) in the same units. 
Could this 'bump' be due to debris produced by the swarm of fragments that was detected 
in the 2006 May close approach? The distribution of infrared-detected fragments from 
Table [3] is shown in Figure [15] directly on top of the 2007 Jan surface brightness profile. 
It is evident that the distribution of fragments (from 2006 May), scaled to the relative 
separation between fragments B and C, is very similar to the debris surface brightness 
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Fig. 15. — Surface brightness profile, at wavelength 24 /im, along the orbit of comet 73P 
in 2007 Jan (solid curve and left axis). The bright fragments B and C are labeled. The 
relative location expected for fragment G is indicated by a dashed line; lack of a detected 
condensation indicates fragment G nearly completely disrupted. The nonzero brightness 
ahead of fragment C is the forward debris trail. For comparison, the histogram (and right 
axis) shows the density of fragments detected in 2006 May. The fragment density is the 
number of fragments with orbital separations, from Table [3X within bins every 17% of the 
separation between fragments B and C. 
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profile; in particular, the high density of fragments, approximately 1.5 times Dbc behind 
fragment C, is very likely associated with the surface brightness 'bump' in 2007 Jan. This 
patch of enhanced brightness is likely a debris field that is much of what remains of the 
fragments that disintegrated during the 2006 perihelion passage. 

It is unlikely that diffuse emission in the 2007 Jan image is due to debris that was 
produced in the 1995 splitting event. This is evident from the sheer amount of material 
required to produce the observed debris field. In the 2006 Apr image, the surface brightness 
along the orbit, between fragments B and C, is approximately 0.2 MJy sr -1 , which 
corresponds to an optical depth of 8.1 x 10~ 10 if we take into account a particle temperature 
of 278/r/j K, where is the distance from the Sun. In the 2007 Jan image, the brightness 
in the same location relative to fragments B and C was 0.9 MJy sr _1 , which corresponds to 
an optical depth of 1.1 x 10 -8 . Thus there is 14 times more material between fragments B 
and C in 2007 Jan than there was in 2006 Apr. (The situation is likely to be even more 
extreme in the region behind fragment B, but this region was not known at the time of 
designing 2006 Apr observations so that portion of the orbit was not observed and cannot 
be directly compared.) 

The mass of debris in the 2007 Jan image is 3 x 10 13 pa cm g, where p is the mass density 
(g cm -3 ) of individual particles and a mm is the particle radius in mm. The particle size is 
likely of order mm size, or else radiation pressure would perturb their orbits significantly; 
however, the 2007 Jan viewing geometry does not completely rule out a contribution from 
smaller particles. The observed width of the debris field perpendicular to the orbit (4.7 x 10 4 
km FWHM), and particle age of 8 months, indicate an expansion velocity of 1.1 m s _1 . If 
the expansion is only due to the momentum imparted to the particles in the gas outflo wing; 



from a nucleus, y e j ~ 2.3(pa cm ) 1//2 m s 1 (jSykes and Walker 



1992 



Lisse et al. 



Reach et al. 



1998 . 



20001 ) . then particle sizes pa cm ~ 4 are inferred for the debris. If the apparent 
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separation of a particle from the projected orbit is due to radiation pressure, a lower limit 
to the particle size can be inferred; we find (3 < 0.01 (or equivalently, pa cm > 0.006). The 
location of the 'bump' in the trail profile also contains clues for the origin of the particles 
and their size. If the particles originated from fragment B in 2006 May, then in order to 
reach their present location at -1.15 Dbc, they must have pa cm ~ 0.006. If the particles 
originate from the large, disrupted fragment G, at original location -1.18 Dbc m 2006 May, 
and they moved to -1.4 D BC by 2007 Jan, then they have pa cm ~ 0.06. If the particles 
originated from the smaller fragment cloud already near -1.4 Dbc i n 2006 May and moved 
less than 0.1 Dbc by 2007 Jan, then we require pa cm > 0.1. 

The observed debris mass of ~ 10 13 g produced by fragment disintegration during 
the 2006 perihelion passage is equivalent to the entire mass of a body 280 m in diameter. 
Fragment B has an estimated 400 m diameter, and fragment G was undoubtedly smaller. 
Thus the mass of the mm to cm-sized debris was comparable to the entire mass of fragment 



G before its disintegration. The debris mass is a lso comparab 



of entire meteoroid streams for other comets (cf. iReach et al. 



e to or larger than masses 



20071 ). However, there are 



clear distinctions between the 73P debris and most other comets. First, the 73P material 
in the 2007 Jan image was primarily produced during the 2006 perihelion passage, whereas 
meteoroid streams for comets with stable orbits are built over multiple revolutions. Second, 
the 73P material originates from small fragments, whereas for most comets the material 
originates from the nucleus. 



8.3. Non-gravitational forces 



Inspecting the available calculations of non-gravitational parameters for other 
periodic comets, we found a trend suggesting a correlation such that comets with high 
non-gravitational forces in their present orbit were more likely to have split in the past. 
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of the JPL So 



ar System Dynamics 



We found this result using the small body database f 

group, together with the historical notes by Kronkf; and Kinoshit43. For all orbits 
with measured non-gravitational forces, we combined the parameters into one number, 
A = 10 9 (Al + A 2 2 + Af) 1/2 . Of 43 comets with A > 3.6, eight had historically observed splits; 
while of the 50 comets with A < 3.6, only 2 had historically recorded splits. For the orbits 
compiled by P. Rochecl, the logarithmic mean A = 17 for 7 split comets, while A = 4.7 
for 65 other comets. Some of the orbits with the highest non-gravitational parameters are 
uncertain, so this result deserves further scrutiny. But the general tendency is unlikely to 
change. 

The elevated non-gravitational forces for fragments of cometary splitting could be 
due to their smaller size or the presence of freshly exposed interior ice, both of which are 
expected properties of fragments. Alternatively, comets with large non-gravitational forces 
may be on the verge of splitting due to greater stresses on the nucleus from the elevated 
outgassing. For 73P, there are orbits with measured non-gravitational terms in 1979, 
1995, and 2001; the values 10 9 {A 1 ,A 2 ) are (6.7,0.4), (6.9,0.5), and (12.8, 1.8), respectively. 
Compared to the ensemble of comets, the non-gravitational forces were moderate in 1979 
and similar during the revolution when the splitting occurred (in 1995). Then in the first 
post-splitting revolution, the non-gravitational forces were much higher. These observations 
support both hypotheses for the correlation between non-gravitational forces and splitting, 
so it is possible that both mechanisms (enhanced outgassing leading to splitting, as well as 
enhanced outgassing from fragments) may be operating. 



3 http: / / ssd.jpl.nasa.gov 
4 http: / / cometography.com 



5 http:/ /www9. ocn.ne.jp/~comet 

6 http: / / www.imcce.fr/page.php?nav=en/ephemerides/donnees/comets/index.php 
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9. Conclusions 

During 2006 March to 2007 Jan, we used the IRAC and MIPS instruments on the 
Spitzer Space Telescope to study the infrared emission from the ensemble of fragments, 
meteoroids, and dust from comet 73P/Schwassmann-Wachmann 3 within a debris field 
more than 3 degree wide. We also investigated contemporaneous ground based and HST 
observations. The Spitzer data spatially resolve the fragments' comae and tails on scales of 
10 3 km. The luminosity distribution of large fragments is similar to the size distribution 
found for nuclei of the ensemble of Jupiter family comets, while there is a paucity of smaller 
fragments. From fragments 73P-B and C, we use the IRAC 4.5 /j,m band to measure the 
fluorescence from outflowing CO2 gas, finding that the abundance of CO2 is 5-10 % of that 
of H 2 0. 

We explain observed morphologies in terms of non-gravitational forces, which we 
quantify using the ratio of solar radiation pressure over gravity, (3, the ratio of momentum 
loss due to sublimation over gravity, a, and the ratio of initial ejection energy to gravity, 
v. Table [6] summarizes the different types of debris and their dynamical influences. The 
large fragments are splayed along the project orbit of the progenitor comet with dynamics 
dominated by v. While the distribution of fragments perpendicular to the orbit is tightly 
peaked to the progenitor's orbit, some fragments are far from the orbit requiring very high 
A% from non-gravitational forces or a extremely large ejection velocity perpendicular to the 
orbit. Very small fragments or very large meteoroids, in the ~ 5 m size range, are seen in 
ground-based and HST images splayed along the dust tail; their dynamics are dominated 
by rocket forces with a > 10~ 3 . Smaller meteoroids, in the mm to cm size range, are seen 
in the Spitzer images only; their dynamics are dominated by weak radiation pressure, j3. If 
the trail particles were icy, they would have large a and would deviate significantly from 
the debris trail. Thus the debris trail particles must be devolatilized. 



Table 6. Summary of fragment and debris size and dynamical influences 



Physical properties Dynamical parameters 

location size ficeZia (3 = F rad / F grav a = F rocket / F grav v = v eo /v orb 



Largest fragments 


along orbit 


200 m 


0.1 


io- 9 


io- 7 


10~ 6 


Fragments 


along orbit 


50 m 


0.1 


io- 8 


IO" 6 


10~ 6 


Mini- fragments 


dust tail 


5 m 


8 


io- 7 


io- 3 


IO" 6 


meteoroids 


along orbit 


5 mm 


< io- 4 


io- 4 


< io~ 4 


io- 4 


dust 


dust tail 


10 /im 


? 


io- 1 


? 


IO" 2 
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